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Abstract. In this work we analyze the validity of Lifshitz’s theory for the case of non-equilibrium scenarios 
from a full quantum dynamical approach. We show that Lifshitz’s framework for the study of the Casimir 
pressure is the result of considering the long-time regime (or steady state) of a well-defined fully quantized 
problem, subjected to initial conditions for the electromagnetic field interacting with real materials. For 
this, we implement the closed time path formalism developed in previous works to study the case of two 
half spaces (modeled as composite environments, consisting in quantum degrees of freedom plus thermal 
baths) interacting with the electromagnetic field. Starting from initial uncorrelated free subsystems, we 
solve the full time evolution, obtaining general expressions for the different contributions to the pressure 
that take part on the transient stage. Using the analytic properties of the retarded Green functions, we 
obtain the long-time limit of these contributions to the total Casimir pressure. We show that, in the steady 
state, only the baths’ contribute, in agreement with the results of previous works, where this was assumed 
without justification. We also study in detail the physics of the initial conditions’ contribution and the 
concept of modified vacuum modes, giving insights about in which situations one would expect a non 
vanishing contribution at the steady state of a non-equilibrium scenario. This would be the case when 
considering finite width slabs instead of half-spaces. 

1 Introduction 

In this paper we are concerned with a first principles cal¬ 
culation of the non-equilibrium Casimir pressure between 
real materials, for the particular geometry of half-spaces 
separated by a vacuum gap of constant width. By non¬ 
equilibrium we refer to a situation in which the objects 
are held at different (but fixed) temperatures. 

Lifshitz formula describes the Casimir pressure in a 
steady situation, in terms of the macroscopic properties 
of the materials. The original derivation [1] was based on 
the use of stochastic Maxwell equations, and relied heav¬ 
ily on the thermal properties of the stochastic sources, as¬ 
suming thermal equilibrium. A few years later, the same 
formula was obtained in Ref. [5] within the framework of 
quantum field theory at finite temperature. Whatever the 
chosen approach, the final result can be expressed in terms 
of the permittivity of the materials (or alternatively in 
terms of their reflection coefficients), and therefore ad¬ 
mits a natural generalization to the case in which the 
half-spaces are maintained at different temperatures, as¬ 
suming that the system reaches a steady state. This kind 
of approach has been followed in Ref. [3], where the non¬ 
equilibrium Casimir force has been considered for the first 


time. Nevertheless, the latter seems rather an extension of 
the original derivation than a full quantum development 
from first principles which is, today, a pending challenge 
in non-equilibrium scenarios. Although there are remark¬ 
able works in the context of Casimir physics (as Ref. 0), 
the present paper focuses on the dynamical aspects of the 
building-up of the steady state of an electromagnetic (EM) 
field interacting with matter in a non-equilibrium scenario 
that, as far as we know, has not been done before, and we 
think that contains relevant physics. 

Beyond the non-equilibrium features, this work also 
enters another (not completely unrelated) issue that has 
been the subject of an intense controversy: whether Lif¬ 
shitz theory is applicable or not to real metals. In the sim¬ 
plest approach, free (conduction) electrons in metals are 
described by the Drude model, although its dissipation¬ 
less limit, the plasma model, seems to be in concordance 
with the experimental results . Some authors claim that 
Lifshitz formula cannot be applied for the Drude model 
[B], while others have reached the opposite conclusion [7], 
claiming that the Lifshitz-Matsubara formula is not cor¬ 
rect for the plasma model. At first glance, our work will 
not address directly this important question, because thro¬ 
ughout the paper we consider a material without free car- 
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riers, i.e., an insulator which permittivity function is given 
by the bounded electrons only (see Ref. [6], for example). 
However, the present work enters the discussion about a 
systematic treatment that allows to understand the limi¬ 
tations of the Lifshitz formula, and how to figure out which 
model is consistent with it. 

A first principles calculation involves a microscopic 
model for the quantum polarization degrees of freedom of 
the materials, interacting both with the electromagnetic 
field and thermal baths that fix the different temperatures 
of the objects (see Ref.[51IH] for more details about phys¬ 
ical aspects of the material model). The evolution of the 
fully quantized system from adequate initial conditions 
should eventually show the emergence of a steady situa¬ 
tion in the long time limit, with Casimir pressure given 
by a Lifshitz-like formula. From a technical point of view, 
this program can be developed in the framework of the 
theory of quantum open systems (see Ref. m)- The quan¬ 
tization procedure can be worked out through Heisenberg 
equations for the quantum operators (canonical formalism 
[TT| 1 or path integral methods (closed time path formalism 
(CTP) [H]), but both are subject to initial conditions m- 
For the latter, after integration of the material degrees of 
freedom and thermal baths, one obtains the influence ac¬ 
tion fRefs. [1^114) 1 for the electromagnetic field. Then, in¬ 
tegrating the field, one can construct the generating func¬ 
tional to calculate the correlation functions (see Ref. m)- 
Using the correlation functions derived, one can compute 
the time evolution of the mean value of the energy momen¬ 
tum tensor, and therefore the Casimir pressure as a func¬ 
tion of time. The main question is therefore whether this 
mean value reaches a well defined limit at long times, and 
whether this limit depends on the initial conditions or not. 
Note that the approaches that are usually implemented in 
Casimir physics can be considered as ‘steady’ quantiza¬ 
tion schemes in the sense that they assume without jus¬ 
tification that the system reaches a steady state and, in 
addition, that the latter is of thermal equilibrium. Both re¬ 
quirements makes the ‘steadiness’ assumption physically 
reasonable and expected to treat equilibrium situations 
(see for example Refs. [T!?l[TH] for the ‘steady canonical 
quantization’ schemes and Ref. m for the ‘steady path 
integral quantization’). Nonetheless, it is not clear how to 
extend it to non-equilibrium (steady) situations or how it 
is done from first principles. As we mentioned before, the 
first work in this direction was Ref. [3] which extends the 
stochastic approach, but a complete quantum approach 
will bring to light the physics of the contributions of each 
part of the model to the electrodynamical quantities. At 
the end, the whole picture gives an appropriate and sys¬ 
tematic way to deal with non-equilibrium situations in 
Casimir physics. 

In previous works, some of us investigated preliminary 
aspects of the problem at hand. In Ref. m the case of a 
quantum scalar field in the presence of an arbitrary mate¬ 
rial was studied. After developing the quantum open sys¬ 
tems approach for the particular situation of a scalar field 
interacting with microscopic degrees of freedom and ther¬ 
mal baths, it was shown that the emergence of a steady 


state independent of the initial conditions is a nontrivial 
issue. While for a quantum field in bulk material such limit 
exists, this is not the case for slabs of material with finite 
width. The generalization to the electromagnetic field has 
been considered in Ref. m, where the CTP generating 
functional for the electromagnetic field interacting with 
a composite environment has been described in detail. 
Specifically, the environment considered there consisted 
of the quantum polarization degrees of freedom at each 
point of space, connected to thermal baths to fix the local 
temperature. Formal expressions for the Maxwell tensor 
and the Poynting vector in terms of the Hadamard prop¬ 
agator were obtained, for the case of the electromagnetic 
field in bulk material, where, once again, the steady state 
is independent of the initial conditions. Some technical is¬ 
sues regarding the gauge fixing procedure in the context 
of CTP have also been discussed. 

In the present paper we will extend the results in [TO] 
to the case of the electromagnetic field in the presence 
of two half-spaces of real media separated by a vacuum 
gap. Considering the previous results mentioned above, 
it is not clear a priori wether the initial conditions will 
be erased or not in the long time limit for this particular 
geometry. Using the formalism of Ref. m we will obtain 
explicit expressions for the Casimir pressure as a func¬ 
tion of time, assuming that the interaction of the vacuum 
field with the materials starts at a given time to- The 
Casimir pressure will receive three different contributions, 
coming from the held’s initial conditions, the polarization 
degrees of freedom, and from the thermal baths. We will 
show that, in the long time limit, only the baths’ contri¬ 
bution survive. This is the contribution expected in the 
Lifshitz approach. Therefore, we will provide a first prin¬ 
ciples derivation of Lifshitz formula for non-equilibrium 
situations in the framework of quantum open systems. 

The paper is organized as follows: In the next Sec¬ 
tion we present the model and the influence action that 
includes noise and dissipation contributions to the calcu¬ 
lation of the Casimir force. Section HI contains the ex¬ 
pression of the Casimir pressure in terms of the electro¬ 
magnetic correlation functions and, following this, the dis¬ 
cussion about the different contributions to the pressure 
appears in Sec. IV. The long-time limit of the complete 
problem is included in Sec. V and finally, in Sec. VI we 
include our final remarks and conclusions. 


2 The model and the field influence action 

In order to include effects of dissipation and noise (fluc¬ 
tuations) in the calculation of the Casimir force between 
two half-spaces of real material interacting with the EM 
field, we will develop a full CTP approach to the problem. 
Therefore, a microscopic physical model for the material of 
the bodies will be introduced. Once the model is defined, 
one expects the material to show effective macroscopic EM 
properties, that will be contained in the Feynman-Vernon 
influence functional included in the CTP formalism. 
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2.1 The Total System: EM Field plus Matter 

As in Ref. [19], we will consider a composite system, con¬ 
sisting in two parts: the EM field (considered as a 
massless vector gauge field) and the real media, which is 
modeled as a continuous set of quantum Brownian parti¬ 
cles representing the polarization density degree of free¬ 
dom P. These degrees of freedom (DOFs) are basically 
three-dimensional harmonic oscillators (with mass 
and frequency 17 [-^1 in each direction) coupled to the field 
at each point of space in a realistic way. The compos¬ 
ite system (EM field and DOFs) is also coupled to ex¬ 
ternal baths of harmonic oscillators in each direction as in 
the well-known quantum Brownian motion theory (QBM). 
Therefore, the total action for the whole system is given 
by: 


Px, q„,x] = + 5o[Px] + ^ ^o[qn,x] (1) 

n 

+ ^Curr[7l'^,Px] +^^i„t[Px,q«.x], 

n 

where the first three terms are the free actions for each 
part of the total system, while the interaction action be¬ 
tween the field and the polarization DOFs is given by 
a current-type interaction which can be written in two 
equivalent ways by defining the conserved current four- 
vector as Jf_i = iy ■ P, —P) and the electric field = 
-djA° - doA^: 


ScurrlA^, Px] = Ao [(<? J^) * = Aq [(<? P^) * , 

( 2 ) 

where A*B = f d^x A(x, t) B{k, t) and Aq is the coupling 
constant between them. It is worth noting that, when writ¬ 
ing the interaction actions, we have introduced the matter 
distribution (/(x), which takes binary values (1 or 0) de¬ 
pending whether or not there is material at the spatial 
point X. 

On the other hand, the interaction between the po¬ 
larization DOFs and each bath is simply linear coupling 
through the constants A^.x for each bath oscillator (n) in 
each point (x) and direction (j): 

S'intpx, q«.x] = J {Pi ■ 9n^x) : (3) 

where A-B = dt A(t) B{t), aJ^Jx are coupling constants 
between the DOFs and the baths, and [j] denotes that this 
superscript is not summed as in Einstein notation. It is 
clear that greek indices sum over 0,1,2, 3 while latin sum 
over spatial components only (1, 2, 3). 

Finally, we will assume that the total system is initially 
uncorrelated, thus the initial density matrix is written as 
a direct product of each part, which we also assume to be 
initially in a thermal equilibrium at a proper characteristic 
temperatures {Pem, Pb,^- the material can also, in 
principle, be thermally inhomogeneous-), 


p{U) = Pem{U) ® ppyi) ® P{q„,,,}(ii)- (4) 


2.2 EM Field Influence Action 

As our goal is the study of the dynamics of the EM field 
correlation functions, we start by coarse-graining the baths 
and the DOFs in order to obtain the exact EM field in¬ 
fluence action. Since we have quadratic actions and we 
choose the initial states as Gaussian, the EM field influ¬ 
ence action will have a quadratic form, with dissipation 
and noise kernels acting over the EM field. 

It is well known from the QBM theory, coarse- graining 
the baths will result in QBM influence actions for each 
component of each DOF in each point of space, i. e., the 
result will be S'if , Px^ ] including the QBM dissipation 
and noise kernels. Therefore, at this point, the EM field 
influence action is defined as: 


J = / dPf / dPi dV[ 


rP'(tt)=Pf 


rP(tt)=Pf 


VV X 


P(ti)=Pi 

/ PP' g'-^o/dx 3 (x)(VA“-P+A-P-VA'“ P'-A' P') 

Jp'(ti)=P[ 

^ g*(So[P]-So[P']-SSip[P,P']) p/^ ^ (- 5 ^ 


Following Ref. m, and considering an initial thermal 
state for the DOFs, the EM field influence action is given 
by: 


Sw[Ai^,Ay = / Pa; / dS' AA>^{x) -2 D^^{x,x') 


X SA^ix') + - N^.{x,x') AA‘'{x') ,(6) 


with AAi^ = A'f^ - A^^, SA)^ = (A'^ -b AA)/2 and the 
dissipation and noise kernels are given, respectively, by: 


D^,{x,x') = r^,^^I5>,k. (7) 

iV^,(x,P) =P^P'= (8) 

which ensures gauge invariance thanks to the fact that the 
differential operator P^i, ° di, ° — 

d^ ^ d, ° +d^ J d, dl, satisfies a'^P^, = a'"P^, = 

0 . 

It is worth noting that the noise kernel presents two 
contributions, one associated to the DOFs and the other 
one associated to the baths: Njk = Nj), -b 

All in all, the dissipation and noise kernels are given 
by: 


lijk{x,x') = Sjk d(x - x') g{x) ^ G|it,x(^ - (9) 
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M-P / i\ X xr /\ f Ppi 

X ) = bjk o(x — X ) -fTi-cotn 


2a 


m 




+ G'iljt - U) g{^ 1 Jt' - U)] , (10) 


Nffc(a;,a;') = 6jk 5(x-x') g{x) • 


r^j] 

'^RetjX 


( 11 ) 

where we are considering materials that can be inhomo¬ 
geneous and anisotropic (birefringent) |19] . 

A^x ^ is the QBM noise kernel generated by the ther¬ 
mal bath over the DOF in the direction j and located 
at X. The Green functions ^ are the retarded Green 
functions for each DOF, which satisfy G[(gj ^^(O) = 0 and 

Gnit x(0) ~ 1-1'^ other words, they are the QBM retarded 
Green functions having x and [_)] as parameters. These 
Green fuctions are defined once the type of baths are cho¬ 
sen at each point of space and direction. It is implicit 
that the chosen directions defining the anisotropic prop¬ 
erties are the Fresnel’s principal axis’ basis (see Ref. HU) 
because, as we shall see, they will allow us to define a diag¬ 
onal permittivity tensor in this basis. In fact, we are con¬ 
sidering the same basis for every point filled with material. 
This is completely general since it is clear, for example, 
that disjoint bodies can have different Fresnel’s basis. To 
include these features, it is necessary to include changes 
of basis between all the present basis. For simplicity, this 
will be omitted and, in fact, for isotropic materials this is 
unnecessary since the permittivity tensor is proportional 
to the identity matrix. 

In the Lifshitz problem there are two different bodies 
interacting, which are assumed to be homogeneous and 
isotropic. Therefore a few simplifications arise. On the one 
hand, homogeneity implies that all the spatial subscripts 
reduce to two labels associated to each of the two bodies. 
Gonsidering two parallel plates separated by a distance I, 
the homogeneity of each body means :x. ^ L,R depending 
on the spatial point lays in the left or right body respec¬ 
tively. On the other hand, isotropy of the material in each 
body implies that superscripts [j] should be omitted be¬ 
cause there is no dependence with the direction. However, 
we will keep generality in the material properties for now. 


3 Casimir Pressure in terms of EM 
correlation functions 


3.1 The Hadamard Propagator and Electrodynamics in 
the Temporal Gauge 

Once we have calculated the EM field influence action 
generated by the material, we proceed to calculate the 
EM field generating functional, which is defined by [T^ : 


r r 

Zc^p[J^.,f^] = dA^ dA>^ dA[^ / VA>^ 

J J JAt-{ti)=A^ 

pA'^(tf)—A^ 

/ ' VA'f^ g^SoM^l-SoM'"]) 

J A'i^{ti)=A'p 

X pEM(Air,Ai;^ti). (12) 

These integrals can be performed using the Eaddeev- 
Popov procedure (adapted to the GTP formalism), to ex¬ 
tract the redundant sums over paths on the same gauge 
class. In m, the gauge condition was introduced in the 
CTP-action as a typical gauge fixing term. Then, the GTP 
integral is worked out by writing the paths as a sum of a 
homogeneous solution Aq{x) (satisfying the initial condi¬ 
tions) plus a shift. The corresponding equations of motion 
for this solution are the ones obtained from the EM GTP 
action (including the EM field influence action), which in 
the case of the temporal gauge reads: 


- df,d^ - ^ t^u'^A''{x) (13) 

-I- 2 J d*x' D^v{x,x') A^i^x) = 0, 

where a is the gauge fixing parameter and is a time¬ 
like four-vector which can be taken as (1,0, 0,0) for the 
temporal gauge. These are four equations for the four com¬ 
ponents of the EM held. 

Therefore, the functional integral can be calculated in 
the Landau gauge (where the gauge hxing parameter a 
goes to 0). All in all, the generating functional reads: 


ZctpISJ, A3] = ( 

\ / i;Ai.i;ni 

where/...\ = f dSAi f dSYii ...HEMiAfAi, 2711^, 

ti], and Wem [A7Ai, ATIi, is the initial EM held Wigner 
functional associated to the EM held initial density ma¬ 
trix pEM(ti) in the temporal gauge (see Ref.[Tn] for its 
dehnition). 

It is worth noting that Aq corresponds to the homo¬ 
geneous solution of the EM equations of motion after im¬ 
posing the temporal gauge, which is written as: 


In the present section we will calculate the EM held GTP- 
generating functional and derive an expression for the 
Hadamard propagator. We shall show how the Gasimir 
pressure can be easily written in terms of this Hadamard 
propagator. 


Ai{x) = - J dx' g^tt(x,x',t-<i) AAf(x') 

-b J dx' ^^g^(x,x',t - <i) A'd7f(x'). (15) 
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^Ret corresponds to the retarded Green tensor result¬ 
ing from Eq. USD after imposing the temporal gauge con¬ 
dition in the Landau gauge, which results in a set of three 
equations of motion plus a residual gauge condition ob¬ 
tained from the equation for (which was erased by the 
gauge condition). The initial conditions for the retarded 
propagator are, as usual. 


x', 0) = 0 , x', 0) = - (5(x - x'). 

(16) 

Once we have calculated the EM field generating func¬ 
tional, we can functionally derive it in order to obtain the 
EM n—point functions and, in particular, the propagators. 

The calculation of the Wightman function is straight¬ 
forward, from which we can read the Hadamard propaga¬ 
tor: 


Qi^{xi,X2) = (A^oixi)A^{x2)) (17) 


+ 


jk 


^/Ret * =1= (f/Ret) (a:i,a^2)- 


This expression holds for every initial state of the field 
and depends on the chosen gauge. Note that the Hadamard 
propagator has two separated contributions. The first term 
is entirely associated to the field’s effective dynamics and 
the initial state. The other contribution is associated to 
the material degrees of freedom represented by the noise 
kernel N, which also splits in two contributions due to the 
composite nature of the material (DOFs plus bath in each 
point of space). 

However, all the dynamics is up to the retarded Green 
tensor in the temporal gauge. As it was shown in Ref. US], 
the choice of the Landau gauge naturally implies that 
A° = 0 in order to avoid divergent terms. Considering the 
temporal component of Eq. m, and imposing the gauge 
condition, one obtains a residual condition over the re¬ 
maining components of the EM field, which can be written 
as: 


V ■ 




(18) 


where the permittivity tensor for the inhomogeneous and 
anisotropic material is given by: 


emr{t-t',X.) = Smr (<5(t - t') + -^O.x ^(x) G|^1 ,x(^ “ • 

(19) 

Since the tensor is diagonal, the expression is given in 
the Fresnel’s principal axes basis. The residual condition 
in Eg. dTSl) is a generalization of the condition considered 
in Ref. [16]. It was shown in Ref. [19] that in the case of 
isotropic and non-dissipative material, Eg. dTSl) reduces to 
the generalized Coulomb condition of Ref. m 

V ■ [e(x) A(x,t)] = 0, 


where the permittivity tensor has been replaced by a single 
function. 

Also in Ref. m it is proved that, in the general case, 
the equations of motion for the spatial components of the 
EM field in the temporal gauge can be written as: 


^2 A 

^ + V X (V X A) + ,, 5(x) A(x, t) (20) 

-I- Ao_x ff(x) dt' ^Ret.x(i - t') ■ A(x, t') = 0, 

where ( GRet.x) = Smk Again, from the fact 

that we could write this diagonal tensor (associated to 
the retarded Green functions), it is clear that the chosen 
basis is the Fresnel’s principal axes basis (in general, the 
tensor would be non-diagonal). It is also remarkable the 
appearance of the third term, which constitutes a finite 
renormalization position-dependent mass term for the EM 
field as the one found in the scalar case in Ref. m- As we 
shall see, this term will be irrelevant in the determination 
of the Green tensor. 

Considering the equation of motion for the EM field, 
the retarded Green tensor Q Ret(x, x',t) can be defined 
as: 


0 = ^ + V X V X ^Ret + Ao_x5(x)^Ret(x,x',t) 

+-^o.xff(x) / dt'^Ret,x(i - i') • ^Ret(x,x',<'), (21) 

Jo 

with the initial conditions given in Eq. (1161) . 

Laplace-transforming the equation, we easily obtain: 

Vx V X^Ret + S^^(s,x)-^Ret(x,x',s) = -I(5(x-x'), 

( 22 ) 

where the Laplace transform of the permittivity tensor of 
Eq.dlSl) is 

■^(S, X) = I -b Aq ,, 5(x) ^Ret.x(s) • 

It is clear that including dissipation it is not possi¬ 
ble to define refractive indexes in the time domain. How¬ 
ever, since the Laplace transform of the permittivity ten¬ 
sor turns out to be diagonal in this basis, the refractive 
indexes can be defined in the Laplace variable’s domain. 
In each direction j, we can define the complex refractive 
indexes 

= 1 -b Ag^^ g(x) ,,(s) 

in such a way that the Fresnel’s ellipsoid is a useful picture 
to describe the material’s anisotropy. 

From Eg. ip^ we can see that the Laplace transform 
of the EM retarded Green tensor turns out to be the 
Feynman’s Green tensor associated to the differential op¬ 
erator V X V X -b ’^(s,x)-, and therefore satisfies 
^Ret(x, x', s) = f/Ret(x', X, s). It is also worth to note that 
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the reality of the EM retarded Green tensor in the time 
domain implies that 


3.2 Transient and Steady Pressure 

We would like to write the pressure in terms of the Hadam- 
ard propagator. In Ref. m, an expression for the expec¬ 
tation values of the components of the Maxwell tensor 
is given in terms of the Hadamard propagator, using the 
point-splitting technique. As we are interested in the cal¬ 
culation of the Casimir force in the Lifshitz problem, a 
simpler expression than the one in Ref. [19] can be achieved. 

Considering the symmetry of the configuration, the 
force between the bodies will be given only by the pressure 
in the perpendicular direction to the surfaces, i. e., along 
the direction parallel to the separation distance /, which 
we will call the z axis. Therefore, in the Lifshitz problem, 
the pressure will be given directly by the zz—component 
of the Maxwell tensor, which can be written for a field 
point xi inside the gap as in Ref. [3]: 


Stt L 


(23) 

where the electric field is given by A® = —OqA^, while 
the magnetic field is S® = (V x A)\ A®-! is the diagonal 
matrix = A^^ = 1 = -A^^. 

Then, using the point-splitting technique and typical 
relations between the different propagators (see [H]), the 
expectation value of the zz—component, can be written 
in terms of the Hadamard propagator: 


PCas (^1) 




Ab 

Stt 


lim 

X2^Xi 


d®" (5^®®® dtA2 




Qtr ixi,x2), 


(24) 


once the propagator is renormalized in the coincidence 
limit. 

It is worth noting that, unlike the one found in Ref.|5] 
which is proposed to correspond to the steady situation, 
the last expression actually is the Casimir pressure which 
emerges at the initial time ti. Therefore, it comprises all 
the transient dynamics of the pressure in the way to reach 
its final value at the steady situation. The pressure can 
depend on time but also on space during the transient 
stage, until finally achieves the steady situation, where its 
value results to be time and space-independent: 


Pcas(a:i) ^ PgL- (25) 

As a final remark, given the splitting of the contribu¬ 
tions to the Hadamard propagator in Eq. (HZI), and also 
that N = -b (depending each one on its tempera¬ 
ture), it is clear that the total pressure can be written in 
terms of three contributions: 


^®Cas(a:i) = Pic{xi) + Pdofs{x 1 , f3p^) + Pb{xi, Pb,^)- 

(26) 

Then, each part of the total system will contribute 
to the Casimir pressure at a given space-time point. The 
main subject of study of the next sections will be to de¬ 
termine which contributions will survive in the steady sit¬ 
uation. 

As we are dealing with an initial conditions problem, 
every time variable is defined in the interval [ti,-|-oo). 
Therefore, we can Laplace transform in each time variable 
ti,t 2 inside the coincidence limit of Ea. (l24l) . Introducing 
a Mellin’s formula for each variable, the second term in 
Ea. d^^ is easily written in terms of the retarded Green 
tensor’s double Laplace transform (note that this term 
only involves spatial derivatives): 


lim dr,di,girixi,X 2 ) = 

X2^Xi 

/*q;i+zoo 7 pa2-\-ioo i 


(27) 


2Tri 


2TTi 


X lim 
X2->Xi L 


irs ilm o o nsm/^ , \ 

; e® (Xi, Si, X2, S 2 ) 


where in the r.h.s. the coincidence limit was taken for the 
time variables and where ai ^2 are taken to define vertical 
lines in the si. 2 —complex planes in such a way that all the 
poles of the integrands taken as functions of si and S 2 are 
at the left of these lines. 

Writing each time variable of the first term of Eq. (IMl) 
in terms of Laplace transforms, we find: 


dtidt2Qir{xi,X2) = 

.ai+*oo 


0:1—200 

Oi +200 


27rz 


ds 


i gSl(tl-ti) 


dt^dt^g^ {xi,X 2 ) 
■s+t2 0H"(xi,si;a;2) 


-200 

- 9t+H®®®(xi,<i;a:2) 

/•ai+zoo 7 /•a2-|-ioo 7 

/ 27ri / 27ri 

«/Oi—200 ^ 02—200 

Si£2(9t+H’"(xi,si;a;2)) - £2(9t+H®®®(xi,ti;a;2)) 
/•ai+ioo 7 7 .a 2 -|-ioo 7 

/ ^ ^ g.i(ti-q) gS2(t2-L) 

Jai—ico 27 rZ Jc22—ioc 27 rZ 

S1S2 0 h’”(xi,si;x 2 ,S 2 ) - Si g^{xi,Si;X2,t2) 

- S2 f?H™(xi,ti;x2,S2)+I/H”"(xi,<i;x2,ti) . (28) 


The second term between brackets does not depend on 
S 2 . Therefore, for that term, the integral over S 2 have the 
integrand which is analytic in all the S 2 —complex 

plane, so the integral over any contour vanishes for t 2 > ti. 
The same happens for the third and last terms between 
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brackets. All in all, we have proved that for the first term 
in Ea. (l24ll : 


hm dt,dt,gS^{xuX2) = 

X2^Xi 


ral+^oo 


_ 1 

27rf 


pa2+ioo 7 
fa2—ioo ^TTZ 


X lim ^h™(xi>si;^ 2 ,S 2 ). 

X2^Xi 


(29) 


Finally, the Casimir pressure can be re-written as: 


P^Cas (yXi ) — 


-1 r 

StT 


q;i +200 


dsi 

2Tri 


pa2-\-ic)0 


dS2 

2TTi 


g(si-|-S2)(tl-ti) 


lim 

X 2 -S-X 1 


0®™(si,S2) f/+'(xi,Si;X2,S2) , 


(30) 


where we have defined the operator 
6'"™(si,S 2) = A*^(siS2 


Last expression for the Casimir pressure seems for¬ 
mally the same to the one found in Ref. [3], however there 
are subtle differences all related to the statement of both 
problems. Eg. (13(11) is not a steady-situation expression (as 
the one in Ref. 0), but it also comprises all the informa¬ 
tion about the transient evolution of the zz—component 
and describes the building up of the Casimir pressure in 
the long time regime. In fact, in Ref. [3], the pressure is 
calculated from the electric field correlation, which results 
to be proportional to a Dirac d-function in the frequen¬ 
cies’ difference, imposed by the steady situation formalism 
(based on the fluctuation-dissipation theorem) in stochas¬ 
tic electrodynamics (SED). Thus, the double integration 
is automatically reduced to one integration by the source 
correlation. In the present work, we calculate the pressure 
from the Hadamard propagator, which is related to the 
EM quantum field correlation and, as we shall see in next 
sections, during the transient evolution it is not necessar¬ 
ily proportional to a Dirac ^-function. Finally, all these 
points are also reflected in the definition of the 0 opera¬ 
tor. In Ref. [3], it only depends on one frequency variable 
which appears as a denominator due to the fact that the 
pressure is calculated from the electric field correlation. 
The present operator is a function of two Laplace vari¬ 
ables which appears as multiplicative factor, since we are 
calculating the full time evolution of the pressure from the 
EM field correlation. 

In the next sections, we will calculate all the contribu¬ 
tions to the pressure and study its time evolution, starting 
from the expression in Eq. o. 


4 Different Contributions to the Casimir 
Pressure 

The formal expression for the Casimir pressure obtained in 
Eg. (1311 . is written in terms of the Hadamard propagator’s 


double Laplace transform. From Eq. (II3 it is clear that it 
has two separated contributions, one associated to the EM 
field’s initial conditions (the first term) and the other one 
associated to the material (the second term). In fact, the 
noise kernel also splits into two parts (N = -|-N^), one 

associated to the DOFs (N^) and the other one associated 
to the baths (N^). Then, the Hadamard propagator’s dou¬ 
ble Laplace transform reads: 


^H^(xi,Si;X2,S2) = (|A()(xi,Si)Ag(x2,S2)^^^ ^^31) 
+ -2i, 2 (f?Ret * (+t'N) * (I/Ret)^) (xi, Si; X 2 , S 2 ). 


Each contribution will be analyzed separately. 


4.1 Initial Conditions’ Contribution 

Provided that the homogeneous solution is given by Eq. 
(HU and the initial conditions for the Green tensor are 
given in Eq. (HU), the Laplace transform of the homoge¬ 
neous solution results: 


+ (x,s) = f dx' {Snlix') - s A+(x')) , 

(32) 

where g^^^{x,x',s) is the Laplace transform of the EM 
retarded Green tensor. 

Therefore, the first term of the double Laplace trans¬ 
form of the Hadamard propagator in Eg. (1311) results: 


+ (xi,Si)Ag(x2,S2)^ 

= I dx! j dx!' e](j'gj(xi,x',si) g’^^{x 2 ,x!', 82 ) 
[aTT/ - siA+j (x') [aTT™ - S2AA™] (x")) 


(33) 


iJAi.sni 


It is clear that calculating the brackets average over 
initial configurations will introduce the EM field’s initial 
state through its Wigner functional. Therefore, this con¬ 
tribution clearly depends on the initial state for the EM 
field. However, instead of proceeding directly to the cal¬ 
culation of the average, it is more convenient to calculate 
the brackets as a quantum expectation value. As the field 
is free at the initial time ti, the quantum expectation val¬ 
ues in the operator formalism and these averages (using 
the homogeneous solution) are closely related. In fact, as 
a particular case of the theory developed in previous sec¬ 
tions, it is easy to show that a free theory (without inter¬ 
actions) verifies: 


(lunl - siA+1 (x') Isnr - (x")\ 

\ L J L J / UAi, 


-Si A[{x');nrix'')-S 2 Arix' 


(34) 
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where the r.h.s. corresponds to the quantum expectation 
value of the anticommutator containing quantum free EM 
field operators in the temporal gauge. In the free EM field 
case, the temporal gauge also implies the Coulomb gauge 
condition (V • A = 0) over the remaining components. 
Therefore, the EM free field operators are: 




dk 


x/2^k(2^)3 


^e^(k,A) 


Ok,A e 


——k-x) 


+ 4,A , 


(35) 


^(x) = 


dk 


.. y/2wk(27r)3 ^ 

— ul pd“kti-k-x) 
“k,A ^ 


^ ia;ke'’ (k. A) 




-k-x) 




where A sums over the transverse electric (TE) and mag¬ 
netic (TM) polarizations, dk,A, ^k a 

photon annihilation and creation operators, (k. A) is the 
j—component of the polarization vectors and Wk = |k|. 

Considering the completeness relation for the polariza¬ 
tion vectors: 


^ e'(k,A)£™(k,A)=d'”^- 


.1 


tk 


(37) 


A=TE,TM 


and an initial thermal state for the EM field we get 


(ak,x ale',A') = “ k') (1 -f N{ujk )), (38) 

(®k.A ak',A'^ = <5 aa' i5(k - k') ^(wk), (39) 

(ak,A ak',A') = (4,A «k',A') = (40) 

where N (wk) = —jy is the photon occupation num¬ 

ber for the initial thermal state of temperature / 3 em- Then, 
after a change of variables, Eq. (EH) can be written as: 


snl - siSA\\ (x') \snr - S2SAf\ (x")') 

11 J / SAuSTli 


coth I X 

2cjk(2^)3 V ) V 2 ' 


X e*‘^4x'-x") (5^52+^2) _ 


(41) 


A(,(xi,Si)Ag(x2,S2) 




r dk 

ul um 

s:lm ^ ^ 

coth 

/^EMO^k 

/ 2wk(27r)3 

^ 2 

1 ^k J 

2 


(42) 

[siS2 + a^k] 


X dx' a:(ilet(xi,x',si) 
x( f dx" t;^™(x 2 ,x",S 2 ) 


Having this result, we can obtain an expression for the 
contribution of the initial conditions to the Casimir pres¬ 
sure. Replacing in Eg. (15(11) by the initial conditions’ 

contribution ( An(xi, si)An(x 2 , S 2 )) it is straight- 

\ /x'Ai.x'n; 

forward that: 


Pic{xi, Pem) = —- 
/3EMWk 


X coth 


r dk 

Ulum' 
dm ^ ^ 

J 2 a;k( 27 r )3 

^ 2 

J 

+*°°dsi r^+*“ds 2 

-200 J Oi2 — 

00 

+ u;^) lim 

X2^Xi 

0 ^"(S 1 ,S 2 ) 


^x' (7;({'gt(xi,x',si) 

dx" e^S(x 2 ,x",S 2 ) e-**' 


(43) 


Remarkably, this expression is quite general. In fact, it 
comprises the full time evolution of the contribution of the 
initial conditions to the Casimir pressure at any time for 
any point of space in a vacuum region, since the bound¬ 
aries appear at the initial time ti. This is the reason why it 
is physically expected that the contribution depends not 
only in time but also in space until reaching the steady 
situation. The information about the specific time evolu¬ 
tion of the problem is encoded in the analytical properties 
of the integrand as a function of si and S 2 , i-e. in the poles 
and the branch cuts present, as we shall see below. 

Moreover, it is important to note that this expression 
is valid for any geometry of the boundaries, including at 
least one vacuum region (where we calculate the pressure). 
The information about the boundaries is contained in the 
Laplace transforms of the retarded Green tensors which 
have to be calculated in a specific situation in order to 
obtain a complete result. In next sections, we will show 
how this works for the Lifshitz problem, deriving the long¬ 
time limit of this contribution. 


4.2 Material’s Contribution 


Therefore, the average over the initial conditions in 
Eg. (1551) results: 


Let us consider the second term in the r.h.s. of Eg. (1511) . 
which is associated to the material. Due to the fact that 
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N = -I- N'®, the contribution splits into two contri¬ 

butions, one associated to the DOFs and the other one 
associated to the baths. However, the first step in Laplace- 
transforming the contribution is the same. As the contri¬ 


bution reads £ 1,2 (0Ret*(9t^t/N)*((yRet)^) (xi,Si;X 2 , 


S 2 ) and the retarded Green tensors depend on the time 
differences, i. e., G^et{x,x') = ,t — t'), then the 

products = 1 = involve convolutions in the time variables be¬ 
tween the noise kernel d^^,N and one of the retarded Green 
tensors. Laplace-transforming is straightforward: 


£1,2 (f/Ret * =1= (CyRet)^) (xi,Si;X 2 ,S 2 ) = 

= y dx J dx.' (y^g^(xi,x, Si) £ 1,2 [9(^4,Nim] (x, si;x',S2) 
X e^™(x2,x',S2). (44) 


A few simplifications arise by considering that the noise 
kernels satisfy 

dl,Nimix,x') = 6im 5{x - x') g{x) 

(see Eas. dTUl) and (ITT|) ). Therefore, taking into account 
the (5-functions: 






coth 


Ppi Dx 


(49) 


nSb(si>S 2) = Ag,, gW „x(si) 

( 50 ) 

The last term in the r.h.s. of Ea. (l46l) is the only one 
associated to the baths. Moreover, it is important to note 
that the second and third terms only depend on one of the 
Laplace variables, while the fourth term does not depend 
on them. 

Is clear that in order to continue the explicit calcula¬ 
tion, one has to define the type of baths which are consid¬ 
ered in each direction and point of space for the specific 
problem. However, we can take a further step in the calcu¬ 
lation without loosing generality. From the QBM theory, 
is clear that for any type of bath, the QBM noise de¬ 
pends on the time differences, i. e., in the time domain 
we have that = Nx\t — t'). Therefore, although 

each time variable is defined in the interval [ti,-|-oo), its 
difference results defined in (— cxd, + 00 ). Assuming that 
the kernel verifies the convergence requirements (this can 
be easily implemented by introducing a cut-off function in 
the QBM noise kernel), it can be written in terms of its 
Fourier transform: 


£1,2 (f?Ret * (^tpN) * (I/Ret)^) (xi,Si;X 2 ,S 2 ) 


= J dx g{x) gd^^{xi,x,si) £ 1,2 

X ^Ret(x2,X,S2), 




(■Sl,S2) 


( 45 ) 


where the presence of g denotes the fact that this contri¬ 
bution comes from the material regions exclusively. 

Gonsidering the definitions in Eqs. (HI and (HU , Lapla¬ 
ce-transforming in both time variables the EM noise kernel 
(taking into account that = 0 = NxB(®i>^i) = 

B(^j> •* 2 ) due to the causality of Gp^^ ^), we obtain: 


£1,2 



(si,S2) = S1S2 n[^|p(si,S 2) - Si N[^|p(si,ti) 
- S 2 NWp(<„S 2 )+N™p(t„<,) 

+ S1S2 n|^|b(si,S2), ( 46 ) 


7 VW(t-t') = ^ fvW(^), ( 51 ) 

where A^^ (<^) contains the dependencies on the tempera¬ 
tures of the baths / 3 x,b- 

Using this general form for the QBM noise kernel, its 
Laplace transform is given by Eq. ([Hn|: 


/ + C)0 7 

00 ^ 


iVx(^) 


(si + 5 w)(s2 — ioj) 


( 52 ) 


Once we have calculated each term of Ea. diHll . the 
Eq. (15(1 gives expressions for each contribution to the Gasi- 
mir pressure in Eg. (1261) . Therefore, we can write: 


where each Laplace transform is given by: 


N™p(si,S2) = 


. (Pp. a 


[j]- 


2G 


[j] 


- coth 


X G£t.x(5l) GL.x(s2), 
=Nl\p{U,Sn) 


Agx^x , ( Ppi Gi 

’ ■ coth —— 


2f2K' 


(^SlS2 + 

(47) 


(48) 

(Sn), 


^Ret,x 


-r 


Bdofs(xi,/3p() = -^ / ^ 




Stt 

X lim 

X2—>-Xi 


27Ti 


' Q!2 —200 


aS2 

27rz 


f A^ 

G-’ (si,S 2 ) J dx g{x)- ’ 


2 m 


X coth 


■/3p(I2W' 

(\ 

2 

U 


( [s?G|^lt.x(si) - 1 ] [s2G11L,x(s2) - 1 


+ S1S2 gL,Js2)) 

X C^Ret(xi,X,Si) g^Px2,X,S2) , 


( 53 ) 
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B Xi, x,B Ja,-ioo 




_2 

27rz 


S1S2 e 




lim 


0 ^''(si,S 2 ) / dx g(x) 


(Si + JW)(S2 — iw) X2-S-X1 L 
X 5Ret(xi,X,Si) G|5et,x(si) G^Ret.Js2) ^^x (w) 

X ^Ret(x2,X,S2)j- (54) 


Note that this expression includes an integral over the 
Fourier frequencies. 

Last results are valid in the most general case, which is 
an inhomogeneous and anisotropic material having local 
initial temperatures at each point. However, these expres¬ 
sions are simplified when considering the interaction be¬ 
tween n different but homogeneous and isotropic material 
bodies of volume F„, at homogeneous temperatures. In 
this case, all the spatial subscripts x (denoting inhomo¬ 
geneity) and superscripts [Z] (denoting anisotropy) have 
to be replaced by a single index n denoting the material 
body. Therefore the integral over the space splits as a sum 
of integrals over each volume Vn , and both expressions can 
be written as: 


Pt 


1 


DOFs 




217„ 


coth 


/ /3Pn ■Gn 


/•ai+ico 7 1-02+100 j 

/ ^ ^^^s^ + S2){t^-U) 

>0.1—100 ^TTZ J oi2—ioo ^TTZ 

( 5 lGRet,n('Sl) 1) ( 52 Gl^et,11(52) l) 
-|- 51^2 11 ( 51 ) Gfiet, 11 ( 52 ) 


X lim 

X2-^Xi 


0^’"{si,S2) [ dx5^g^(xi,x,si) 
JV„. 

X ^Ret(x2,X,S2) , (55) 


-I n -\-00 7 p' 

Pb(xi,/3„.b) = 

f02+^00 ^^s^+S2Kt^-t 


o^+^oo 


_ 1 

27ri 


27Ti (si + iu;)(s2 — iuj) 


11 ( 51 ) Gfiet, 11 ( 52 ) 


X lim 

X2—^Xi L 


0^^isi,S2) I dx t;]^gj(xi,x. Si) 

JV„ 


X eRet(x2,X,S2) . 

In this case, the quantity 


(56) 


lim 
X2^Xi L 


0^'"{S1,S2) I dxg;^'^t(xi,X,Si) C?^'3t(x2,X,S2) 


'Hi 


is common to both contributions. However, the Mellin’s 
integrals are the ones that define the time evolution and 


steady state of these contributions, and are very different 
due to the analytical properties of each integrand. 

Moreover, is important to have in mind that in these 
expressions, it is implicit that at every point of the ma¬ 
terial bodies, the basis of Fresnel principal axis is the 
same. Each point defines a Fresnel’s ellipsoid having its 
axis along the directions in which each component of the 
DOFs fluctuates. This is, in principle, a limitation in the 
present model, although it could be easily fixed by consid¬ 
ering different basis in each body (or even at each point of 
space). As we mentioned before, in order to simplify the 
calculations but without losing the anisotropy properties, 
we will keep an unique basis for all the material points. 
This limitation disappears when considering isotropic ma¬ 
terials, as in the Lifshitz problem analyzed below. 

However, it is important to recognize that the ellip¬ 
soid’s picture is valid for real refractive indexes. In our 
case, as we shall see, due to the inclusion of dissipation 
in the problem, the three refractive indexes are complex, 
albeit the Fresnel’s ellipsoid results useful as a conceptual 
and practical picture. 

As final remarks, on the one hand, it is worth noting 
that considering this case for the initial conditions’ contri¬ 
bution has no formal simplification over Eq. (1431) since the 
initial field modes and the spatial integrals are over the 
whole space (^(x) does not appear in those integrals). On 
the other hand, it is important to take into account that 
all the expressions containing a coincidence limit have to 
be regularized before calculating the limit. For the case 
of the Lifshitz problem, this will be done after integrating 
over the space. 


5 The Long-Time Limit of the Full Problem 

With the general expressions for each contribution of the 
Casimir pressure generated by different bodies of homoge¬ 
neous and isotropic material at different temperatures in 
a vacuum region, we can study the Lifshitz problem con¬ 
sisting in two half-spaces by letting n be L or R for the 
left and right plate respectively, as we have mentioned at 
the end of Sect. 12.21 


5.1 Analytical Properties of the Green Functions for 
the Lifshitz Problem 

From Eas. (H51) . (1551) and (1551) the time evolution and steady 
state of each contribution will be governed by the analyti¬ 
cal properties of the integrands as functions of the Laplace 
variables si, S 2 - Therefore, for a given problem is critical to 
know the analytical properties of both Laplace transforms 
of the QBM’s retarded Green function GRet,n(s) and EM 
retarded Green tensor s). 

As the DOEs are considered as quantum Brownian 
particles, irrespective of what type of environment is con¬ 
sidered the Laplace transform of the retarded Green func¬ 
tion can be easily obtained from the equation of motion 
of the QBM’s theory. In the present case, the result is the 
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same as the one found in Ref. m- Then, we can easily 
write: 


GRet,„(s) = 7 - 7 , (57) 

+ 122 _ 2 

where Dn{s) is the Laplace transform of the QBM’s re¬ 
tarded Green function for the environment (baths) of the 
plate n. 

Thus, given a spectral density for the baths, the lo¬ 
cation of the poles will define the time evolution and the 
asymptotic behavior of the retarded Green function. How¬ 
ever, causality implies, by Gauchy’s theorem, that the 
poles of GRet,n should be located in the left-half of the 
complex s—plane. In fact, assuming that 12„ ^ 0 and that 
the baths include cutoff functions in frequencies, the real 
parts of the poles are negative (see Ref. [15] 1. 

On the other hand, the EM retarded Green tensor is 
defined by the equations of motion obtained from the EM 
CTP action of Ea. (IT51) after imposing the temporal gauge. 

In the inhomogeneous and anisotropic case, by solving 
Eo. d^^ for a given permittivity tensor defining boundaries 
and material bodies, we can obtain the Laplace transform 
of the EM retarded Green tensor and calculate the Casimir 
pressure through Eos. d^ . (l53ll and (l54)) . However, the 
solution in these cases is very complicated. 

Therefore, we turn to the simpler case of considering 
the Lifshitz problem, where two different homogeneous 
and isotropic parallel half-spaces are separated by a vac¬ 
uum gap of length I along the z direction. The origin of the 
coordinate system is defined in the middle of the gap. In 
this case, each contribution to the Casimir pressure has to 
be evaluated from Eqs. (021), (1221) and (1551) . where the last 
two equations are simplified versions of the general inho¬ 
mogeneous and anisotropic case, but the first one remains 
unaltered. In this case, the Fresnel’s ellipsoid turns out to 
be a sphere and the permittivity tensor results a function 
times the identity matrix, allowing any basis to describe 
each material body. Therefore, as we anticipated before, 
we omit the superscripts [j] and replace x by the indices 
L, R related to each material body. All in all, Eo. d^^ sim¬ 
plifies to: 

^VxVx-b e(s, z)^^Ret(x,x',s) = -I5(x-x'), (58) 

where the refraction index only depends on z and it is 
given by 

n2(s, z) = e(s, z) = 1 -b ® G'^et(s) + 

^0,R ® ~ 2 ) ^Ret('S)- 

It is remarkable that this equation and the one solved 
in Ref.|3] are basically the same, but this time the solu¬ 
tion gives the Laplace transform of the retarded EM Green 
tensor. Two differences are in order: the equation is for¬ 
mally the same by replacing the Fourier variable in the 


solutions given in Ref. [3] by is; on the other hand, after 
the replacement, the r.h.s. of the Green function equation 
in Ref.jSj is the r.h.s. of Eg. (1551) times dvrs^, because it 
is compensating the differences in the definitions of the 
operator G®"* in the pressure expression of Eg. (ISO)) as it 
was commented at the end of Sect l3.2l As the Laplace 
variable appears as a parameter in the equation, we can 
easily obtain the Laplace transform of the EM retarded 
Green tensor by dividing the one given in Ref. [3] by 47rs^. 
Due to translational invariance in the parallel coordinates, 
the solution is given in terms of the Fourier transform in 
those coordinates: 


(59) 

where Q = {Qx,Qy, 0 ), being x,y the parallel directions. 


For a held point inside the gap (—| < z < ^), the 
transform of the EM Green tensor depends on the value 
of the source point z'. For z' < — | we have: 


A 
D,, L 




2(Zz ^=TE,TM 

+ 6/5,4-] e-Me^Y+] (60) 


and for z' > A 


£?44t(2^:2;',Q,s) = 


E 


+1^ r 


-1 

2 „( 2 ) ^ 

zyz ^=TE,TM ^ 




+ 6^,4+] (61) 


while for a point source inside the gap — ^ < z' < | the 
tensor splits in bulk and scattered contributions = 

Qlj I Qij \ . 

^Ret,Bu ' ^Ret,Sc/' 


^Ret,Bu 

1 


(z,z', Q,s) = - 


Siz Si 


S{z — z') 


2qz 




-q^iz-z') Ci 


0{z 


/j=TE,TM 


6/1,51 I 6u,il I e 




0(z' - z) 


(62) 


^Ret,Sc(2^> 2 ^ > '®) — 9 E/ n 6 / 1 , 4 +] 6 ;ij-[+] 

//=TE,TM ^ 

X e-«A--'+20 ^y^.^^iz+z'-i) 

+ e,A-] e,A+] +e,.4-] e,.,]-] 


X e 


qz{z — z' —21) 


(63) 


Here we are considering the notation given in Ref. [3] 
adapted to our case. Therefore, the EM (complex) wave 
vector in the medium n is given by: 
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q(-)[±]=Q±*9i")z, (64) 

with n = L,R corresponding to each plate (while omit¬ 
ting the indices for the vacuum region) and where the sign 
-I- (—) corresponds to an upward (downward) wave. The 
vector Q is always a real vector and appears as the pro¬ 
jection of the wave vector [±] on the interface, while 
the 2 ;—component is given by: 

= \/en{s) + < 52 . (65) 

The wave vector q^"^ [±] lies in the plane of incidence 
defined by Q and z. We can introduce the tranverse elec¬ 
tric (TE) and magnetic (TM) polarization vectors: 

^teW = Q X z, (66) 


(n) 

®TM 


[±] = e 


(n) 

TE 


[±] X q(”) 


[±] 


Q iT Q 

y/snis) is 


(67) 


The reflection (r) and transmission (t) Fresnel coeffi¬ 
cients for each single surface and components are given 
by: 


r 


TE 

n 


qz 


(n) 


qz 


(n) 

qz 


TM 
' n 


(n) 

gn Qz-Qz 

I ('^) ’ 

£n qz + qz 


( 68 ) 


t 


TE 

n 


2 gi 


(n) 


I (^) 

Qz + qz ^ 


^TM 


2y/£n{s) gi"^ 
£n qz + gi"^ 


(69) 


while the multiple reflections (which does not enter in the 
bulk part) are described by the denominator: 


D^ = l-r^ ri^ (70) 

Beyond the apparent complicated expressions for the 
Laplace transform of the EM retarded Green tensor, its 
analytical properties can be analyzed relatively easy. Each 
expression presents several poles and branch cuts. The 
poles basically are the ones associated to and there is 
one pole located at the origin s = 0 . 

The poles coming from should be calculated in or¬ 
der to obtain the time evolution of the system. However, 
causality implies that all these poles must lie on the left- 
half of the complex s—plane in order to give convergent 
contributions in the evolution. It may happen that some 
poles have vanishing real part. These poles would con¬ 
tribute to the long-time behaviour of the EM retarded 
Green tensor (and therefore to the pressure). However, if 
we set s = i Im(s) (with Im(s) ^ 0 ) it can be easily shown 
that there is no poles of this type. 

The pole at the origin appears explicitly in the first 
term of Ea. (l62ll . However, in the second term and in Eqs. 
dSOl), (El) and the pole is only associated to the terms 
with fj, = TM, since in any region the TM polarization 


vectors provide the pole through its denominators, as is 
clear in Eq. (El). 

Besides the poles, the Laplace transforms in each re¬ 
gion present several branch cuts, which provide to the 
transient time evolution of the EM Green tensor. In the 
present case, by inspection, the branch cuts are given by 
the square roots + /c| and ^en{s)s‘^ + fc|. The extra 

branch cut e„(s) that appears in the expression of the 
transmission coefficient for the TM terms in Ea. dS^ that 
enter the Laplace transforms in Eqs. dSOl) and (El) should 
not be considered, since it cancels with the same square 
root contained in the denominator of the polarization vec¬ 
tors e^(±) of Ea. dHT)) . that also enter the same Laplace 
transforms. 

Summarizing, the analytical properties of the Laplace 
transform of the EM retarded Green tensor are very sim¬ 
ple, despite of its complicated expressions. These proper¬ 
ties are important in order to compute the full time evo¬ 
lution of the problem and determining the long-time limit 
of the different physical quantities of interest. 


5.2 The Steady Situation of the Lifshitz Problem 

Eqs. (El-dSSl) contain the Laplace-Fourier transforms for 
the EM retarded Green tensor for a field point inside the 
gap in the Lifshitz problem, from which we can calculate 
the corresponding Casimir pressure in the steady state. 

For all the three contributions of Eos. (HU . (IHK1) and 
(ED, in order to completely solve its time evolutions we 
have to calculate both inverse-Laplace transforms through 
the residue theorem and the analytical properties associ¬ 
ated to branch cuts. 

However, as our aim is to determine the long-time limit 
of the contributions, the work is simpler and related to a 
subtle issue of the long-time regime. The ‘steadiness’ of the 
physical quantities as the pressure or energy is related to 
the fact that they do not depend on time in the long-time 
regime. In other words, usually for the systems of interest, 
due to physical/intuitive reasons, the complete system is 
expected (or assumed) to be in a steady situation in the 
long-time limit. In many cases, this allows one to face and 
study the steady regime directly, without a rigorous (and 
unnecessary) derivation from a full dynamical problem. 
In the present case, we are in the opposite situation. We 
have solved the full dynamical problem in terms of double 
inverse-Laplace transform, that should be calculated using 
the residue theorem. But we cannot do that completely 
due to the technical impossibility of locating the poles of 
the integrands in the complex planes analytically (we only 
know a few properties related to causality of the retarded 
Green tensors). 

Nevertheless, to obtain the long-time limit of quanti¬ 
ties given in terms of double inverse-Laplace transforms 
depending on jg relatively easier. 

As the time evolution depends on the configuration 
of poles and branch cuts that the integrand presents, we 
have to distinguish which analytical properties contribute 
to the steady state of the quantities. As it is stressed in 
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Ref. [20] ; we can show that the typical branch cuts that ap¬ 
pear in the Laplace transform of the EM retarded Green 
tensor do not contribute to the long-time regime. Then, 
the steady situation must be given by some of the poles, 
which always result in exponential temporal behaviours 
depending on the poles’ values. For example, consider¬ 
ing a typical Laplace transform (i.e., depending only on 
one Laplace variable), if a given pole have a non-zero real 
part, through causality we know that it must be a negative 
number. Therefore, the resulting behaviour associated to 
that pole will be an exponential decaying with time. On 
the other hand, if a given pole has vanishing real part (i. 
e., the pole is purely imaginary), the resulting behaviour 
will be oscillating in time with a frequency given by the 
imaginary part of the pole. As a last case, it turns out 
that if the pole is 0, the resulting behaviour will be of a 
time-independent constant. 


All in all, in order to determine the steady state of 
each contribution it is necessary to study the combination 
of the poles present in both integrands over si and S 2 - 


5.2.1 The Long-Time Limit of the DOFs' Contribution 


We can start by considering one of the parts of the mate¬ 
rial’s contribution to the Casimir pressure for the Lifshitz 
problem, which is associated to the DOFs’ contribution 
given by the Ea. dHKI) where, as we anticipated, n = L,R. 
The quantity 


lim 

X 2 ->Xi 


0^''{si,S2) [ dx^^^^t(xi,x,si)0^tt(x2,x,S2)], 
Jv„ J 


However, in our case, the presence of 
in the double inverse-Laplace transforms makes the poles 
associated to si to combine with those associated to S 2 
in order to determine the time evolution of the quantity. 
Thus, the case of a pole at s = 0 is not the only way 
we can obtain time-independent behaviours in our quan¬ 
tities of interest. In other words, if we solve both Laplace 
integrals, each term of the full time evolution of the quan¬ 
tity will be a combination of two of the poles, one associ¬ 
ated to the Laplace integral over si (for example Si S C) 
and the other one to S 2 (consider it at S 2 = S 2 )- Then, if 
we combine two poles satisfying Si = —S 2 , there will be 
no exponential associated to that term in the time evo¬ 
lution, obtaining a time-independent term. In particular, 
is clear that the case of both poles equal to 0 satisfies 
the required condition, giving indeed a time-independent 
term. However, if we take a pole with negative real part, 
for example Si in such a way that Re[si] < 0, the re¬ 
quired condition to obtain a time-independent term would 
S 2 = —Si. As this implies Re[s 2 ] > 0, which is forbidden 
by the causality property, then for all physical quantities 
builded from causal quantities, all the poles presenting 
negative real part cannot be combined in order to obtain 
a time-independent term. 


which is present in this contribution as well as in the baths’ 
one can be simplified in the present case, since in the last 
Section we have given the Laplace transform of the EM 
retarded Green tensor for the Lifshitz problem and its 
analytical properties. 

Taking into account that for this problem the parallel 
coordinates can be Fourier-transformed as in Eq. (1551) . the 
integrations over the parallel coordinates xy can be easily 
done, obtaining: 

/ dyi C?](.gj(xi,x,si) eRet(x2,x,S2) = 

JVr, 

= / (^ dzgl,{z„z,Q,s,) 

X Ct(^2,^,-Q,S2), (71) 

where in the r.h.s. the integration over Vn implies the in¬ 
tegration over each half-space. Then, we can write 

r /■(-i)"o<3 

/ dz = (—1)" / dz, 

Jv„ 


As a last possibility, let us consider the poles with real 
part equal to 0. In these cases, considering a pole isi with 
Si G R, we can verify the required condition with a pole 
iS 2 = —isi in such a way that both oscillating evolutions 
cancel out in the exponential giving a time-independent 
term. 


where we associate n = L [n = R) on the l.h.s. with n = 1 
(n = 2) on the r.h.s. of this equality. 

By considering Eqs. deoj) and (EH, is clear that the re¬ 
maining integration over z involves exponential functions 
in the second argument, being easily calculated, obtaining: 


Following this train of thought, the time-independent 
terms will be so from the very beginning. The steady sit¬ 
uation rises when the transient term vanish, i. e., when all 
the time-dependent terms resulting from poles combina¬ 
tion and branch cuts go to 0. Thus, the relaxation time of 
the system will be defined by the last term that vanishes 
between all the transient terms. Moreover, the relaxation 
time will be equal to the smallest non-zero real part of 
the combination of poles or branch points in the three 
contributions. 


/ rfx C;^^gt(xi,x,si) C?^tt(x2,x,S2) = (72) 



X ^^tt(^i>(-irv2,Q,si) gitAz2,{-m/2,-Q,s2). 


Therefore, from Eq. (1551) , we finally obtain the full time 
evolution of the contribution to the Casimir pressure: 
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Pdofs(xi,/3p„,/) = -^ ^ 21? 

n=L,R 


\,n^n 


coth 




dQ 


na.i-\-ioo 


dsi 

2Tri 


^ 012+200 


ds 


£ p{si+S2){tl-ti 


2TTi 


ev-.-.v- (73) 


(Sl GR,et,n(si) ~ l) (^2 GRet,n(s2) — l) + ‘!^^SlS2 


G^Ret,n (^1 )G^Ret,n (^2 ) 


lim 

X 2 ->Xi 




i-m/2, Q, s,) i-iri/2, -Q, 52) 

(q'Y\si,Q) + qY'\s2,Q)^ 


where we have stressed the fact that the contribution to 
the pressure depends on the plates separation 1. 

Now, in order to determine the long-time behaviour 
of this expression we have to analyze the integrand’s an¬ 
alytical properties in terms of Si and S 2 simultaneously. 
For a given n, all the terms will not contribute necessarily. 
Taking in account the analytical properties of the retarded 

Green functions Gnet.n and the EM retarded Green tensor 
^^ 

G Ret (which splits into TE and TM terms) commented 
in the last Section, by inspection it turns out that the only 
poles of the whole two-Laplace variables integrand result¬ 
ing in a time-independent term are the ones located at the 
origin simultaneously for both variables. In other words, 
the poles located at the origin for both Laplace variables 
(si = 0 = S 2 ) associated to the products TM x TM re¬ 
sulting from 

i-mn, Q, si) Ct(^ 2 , {-iri/2, -Q, S 2 ). 

(74) 

However this is not so straightforward because this prod¬ 
uct is also multiplied by 

[(SlGRet,n(si) — l)(s2^Ret,n(s2) — 1) + 


■G^SiS 2 GRet,n(si)GR, et,n (S2)] 

and 0J'=(si , 52 ), containing linear factors si and S 2 that 
can eventually prevent of si = 0 or S 2 = 0 of being poles 
in a given term. At the end, we find poles of first and 
second order depending of which terms are considered. 

As we mentioned before, the poles at the origin are 
related to the terms associated to the TM polarization 
vectors of Eq. (EZD- Therefore, from Eqs. (EOl) and (| 5 T|) , we 
can see that the term TM x TM of the product in Eq. (TMl) 
contains 1/sf and l/s| as factors. 

For example, terms associated to the combinations sf 
GRet.n(si) (or s| Gnet, 11 (^ 2 )) will cancel out the denomi¬ 
nators 1/si (or I/S 2 ) and therefore those terms will have 
no pole at Si = 0 (s 2 = 0). Those terms contribute only 
to the transient regime. 

On the other hand, the term independent of si and S 2 
resulting from 

(s? GRet,n(si) ~ 1)('S2 GRet,n(s2) — 1) + 


I 7 yjSiS 2 Gflet,n(' 5 l)Gpiet,n(' 52 ) 

will give poles of first order for both Laplace variables 
through it combination with the first term of 0 -^^(si,S 2 ) 
(associated to the electric field’s contribution), and poles 
of second order for both variables in the combination with 
the second term of 0^^(si, S 2 ) (associated to the magnetic 
field’s contribution). 

Finally, the term G^siS 2 GRet,n(si) GRet.n(s 2 ) will only 
present poles of first order when combined with 0^^ (si, S 2 ) 
for the product with the second term (the combination 
with the first term will not have poles at the origin). 

Therefore, the contribution to the pressure in Eg. (1751) 
can be rewritten by grouping the terms according to the 
order of the pole at the origin in both variables: 


Pdofs{xi, = (Terms with second order poles at 0) 

-I- (Terms with first order poles at 0) 

-I- (Terms without poles at 0). (75) 

In first place, it is clear that the last term has no con¬ 
tribution to the steady situation, taking part only in the 
transient stage. 

In second place, it is easy to show that the second or¬ 
der poles at the origin result in terms that are directly 
zero or that diverge in the long-time limit (ti —^ — 00 ). 
This is because, in the residue calculation of the second 
order poles, it is necessary to differentiate the integrand 
with respect to the Laplace variable in consideration (ei¬ 
ther Si or S 2 ). At this point, it is important to consider 
that the residue in these terms can be analyzed separately 
in each variable and that the integrands have the form of 
gs(t-ti) ^jg| 2 es a function which depend on \/s^+C^ and 
\/£„(s) s^ -I- Q^. First, differentiating the exponential re¬ 
sults in times the same function depending 

on y/s^ + Q^ and y^£„(s) s^ -I- Q"^. The evaluation of this 
on s = 0 gives the mentioned term that diverge in the long¬ 
time limit. However, all these terms always have another 
term which cancels them, giving no contribution. Sec¬ 
ondly, differentiation the function depending on ys^-l-t^ 
and yj£n{s) gives terms that appear accompanied 

by s/\/s 2 + Q 2 0 ,- [£(^(s) 5 2£„(s)]s/2y'£„(s) s^ -t Q^, 

which evaluated on s = 0 are the mentioned vanishing 
terms. 

In conclusion, the second order poles at the origin do 
not contribute to the steady situation. 

The remaining terms are the ones containing first order 
poles at the origin. These terms, in the temporal and spa¬ 
tial coordinates’ domain, are associated to a time deriva¬ 
tive of one of the Heaviside functions 0(t — ti) which 
appear in the definition of the retarded Green functions 
or the EM retarded Green tensor due to its causal be¬ 
haviour. Mathematically, deriving a Heaviside function 
gives a Dirac J-function. Therefore the expression is pro¬ 
portional to ^(0). However, this is caused by the sudden 
beginning of the interaction. Physically, this is an approx¬ 
imation to the fact that the interaction rises in a very 
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short time. If instead, for example, we replace the Heavi¬ 
side function by a smooth function going from 0 to 1 in fi¬ 
nite time, approaching this value asymptotically, we would 
avoid this problem, keeping the convergent behaviour of 
the whole expression at the long-time limit. In summary, 
the first order poles at 0 are the result of deriving a Heav¬ 
iside function representing the sudden switching-on of the 
interaction and hence they are not physical. Therefore, 
they will not have contribution to the long-time regime. 
This agrees with the limit for an analog contribution con¬ 
sidered in Refs. mm, where the same contributions in 
the coordinates’ domain are discarded by the same rea¬ 
sons. 

Then, we have proved that for the contribution of the 
DOFs to the pressure, the long-time limit {ti —>■ —oo) is 
given by: 

Fdofs(cci,/3p„,0 ^ F^ofs(/3p„,0 = 0- (76) 


5.2.2 The Long-Time Limit of the Baths’ Contribution 

Now we proceed to calculate the contribution resulting 
from the baths. 

We recall the general expression of Eg. (1551) found for 
the baths’ contribution to the pressure in a scenario of 
homogeneous and isotropic bodies, considering as in the 
last case n = L, R. 

As in the DOFs’ contribution, the pressure contains 
the expression 

/ C?^^gt(x2,X,S2)- 

For the case of half-infinite parallel plates, Eq. ([7^ pro¬ 
vides a simplification for the expression and then, the pres¬ 
sure reads: 


-| _ p-\-00 7 

Fb(xi,/3„,b,0 = -^ 


rvi — 7 ('vi StTZ , 


n=L,R 

■^ 2 -ioo 27rj (si -h iuj){s 2 - iuj) 


(qi''\si,Q) + qi"\s2,Q)^ 


lim 

X2^Xi L 


e^Hsi,s2) 


X gltetiz2,{-iri/2,-Q,S2) . 


dQ 

(27r)2 


(77) 


Despite the branch cuts, in order to determine the 
long-time behaviour of this contribution we have to con¬ 
sider again combinations of poles in the Laplace variables 
Si and S 2 in such a way to obtain a time-independent 
term. 

As a first step, as in the DOFs’ contribution, we have 
to analyze the pole at si = 0 = S 2 coming from the TM x 
TM terms. However, due to the presence of a factor si S 2 


in the integrand, the only term of 0 ^^(si,S 2 ) that can 
have a contribution through this pole is the second one 
y^pm ^prj^mik (associated to the magnetic fields), 

giving first order poles for both Laplace variables. On the 
other hand, the first term of 0^^(si, S 2 ) (associated to the 
electric fields) contains an extra factor si S 2 which finally 
cancels the denominators in the terms TM x TM which 
provide the pole at si = 0 = S 2 , having no pole for the 
first term. 

The calculation of the pole implies summing the in¬ 
dexes in the expressions. Calculating explicitly all the sums 
and taking into account the coincidence limit, it is easy to 
see that the numerator also vanish for si = 0 = S 2 . Then, 
when the denominator vanishes, the numerator vanishes 
too. Using L’Hopital rule one obtains a finite value, and 
therefore, si = 0 = S 2 it is not a pole. 

Besides the point si = 0 = S 2 , by inspection is clear 
that this contribution also presents another combination 
of poles on si and S 2 that produces a time-independent 
term. Due to the fact that the QBM noise kernel depends 
on the time difference (and then it can be written in terms 
of its Fourier transform), the denominators (si + *w) and 
(s 2 — iuj) appear in this contribution denoting the non- 
damped dynamics of the baths. Therefore, the pole at 
Si = —iuj combined with the pole S 2 = iui gives a time- 
independent term. 

All in all, the long-time contribution of the baths to 
the Casimir pressure can be written as: 


^ + 00 


doj 

2tt 


UB“(/3n,B,0 = ^ E / 

^ GRet,n( 

(-iw, Q) + gi"' (iw, Q)^ 

X e*Q-(-qi-^ii) e^'^^(zi,(-I)";/2,Q,-*a;) 

X 0^tt(^2,(-l)”V2,-Qua;) 


Uj‘^Nn{uj) 


dQ 

w? 


lim 

X2-i-Xi L 


ioj) 


(78) 


This is the contribution of the baths to the Casimir 
pressure in a general non-equilibrium context. An impor¬ 
tant connection with previous works can be established by 
employing a fluctuation-dissipation-type relation into the 
result. Besides that it is well-known that the fluctuation- 
dissipation theorem is valid for systems in thermal equi¬ 
librium, and the total system in this case is out of equi¬ 
librium, we can use a fluctuation-dissipation-type relation 
for the dissipation and noise kernel of the QBM (D„ and 
Nn respectively), i. e., the kernels generated by the baths 
in each point of space acting over each DOF. Therefore, 
we can write, for the Fourier transforms of the kernels, the 
fluctuation-dissipation relation: 


Nn{uj) = coth Im [D„(a;)] , (79) 

where Dn{uj) is the Fourier transform of the QBM’s dis¬ 
sipation kernel. 

From the definition of the permittivity tensor in Eq. 
(HU), it can be proved that for the present case, the Fourier 
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transform of the permittivity function is given by 
£ra(w) = 1 + Aq „ GRet,n(w), 

having that G'Ret,n(—w) = Gj^et n(‘^) reality of 

the QBM Green function. On the other hand, since the 
Laplace transform of the QBM GRet,n(s) of Eq. (EB is as¬ 
sumed, by causality, to have poles with n^ative real parts, 
then it is verified that GRet,n(—*w) = GRet(w) and the 
same happens for the QBM dissipation kernel Dn{—iui) = 
Dn{ui) (in fact, the connection between the Laplace and 
Fourier transforms applies for every causal function which 
is 0 for null value of its variable). Therefore, it is straight¬ 
forward to prove that: 


in Ref.|3] for the steady situation of the Lifshitz prob¬ 
lem. However, following this procedure, it is worth noting 
that this result can be extended for the case of inhomoge¬ 
neous and anisotropic materials since both purely imagi¬ 
nary poles are always present as it can be seen from the 
general expression of Ea. dMl) provided one can calculate 
the Laplace transform of the EM retarded Green tensor 
for a given problem, which is the main difficulty in most 
cases. Then, we can in general write: 




X ^Ret(xi,x,w)c0th 


2tt 
/3x,b 


©•’'"(a;) / dxg(x) 




Im 


w 


Im[£„(a;)] = Im |73„(a;)] |GRet,n(w)|^ (80) 

= Im [D„(a;)] GRet,n(-*w) GRet,n(*w). 

Then, introducing Eq. Ci) into Ea. dTRll and using Eq. 
dlol), we obtain: 



(83) 


which is the full generalization of the steady result found 
in Ref. [3]. 


P^{Pu,b 3) = -^ Y, f ^w'coth 


X Im [e„(w)] 


Stt 

dQ 


n=L.R 


/3n.BW 


X lim 

X 2 — >-xi L 


(27’')^ ^q'f-\-iuj,Q) + 


(81) 


(-1)"V2, Q, (-1)”V2, -Q, 


At this point, we use Ea. d73|) in order to go back to 
the spatial integral having si = —ioj and S 2 = ioj. Then, 
we can firstly employ the property associated to the real¬ 
ity of the EM retarded Green tensor in the time domain 
for the last factor t7^g^(x,x', zo;) = , —ioj), fol¬ 

lowed by the property of being a Feynman propagator to 
the last Laplace transform of the EM retarded Green ten- 
sor t/^(fj(x, x', —zw) = 0^^g).(x', X, —zw). Finally, for both 
Laplace transforms of the EM retarded Green tensor we 
use the connection between the Laplace and Fourier trans¬ 
form ensured by the causality behaviour to obtain: 





xlm[e„(a;)] lim dx 

X 2 ->X 1 L Jy^ 

X eRrt(x,X2,a;) 


(82) 


which is exactly the result obtained in Ref. [3] as the to¬ 
tal pressure for the Lifshitz problem, where we have set 
= 0^'^{—iuj,iuj) from its definition. 

All in all, we have proved that the baths’ contribu¬ 
tion in the long-time regime gives exactly the result found 


5.2.3 The Long-Time Limit of the Initial Conditions' 
Contribution 

After determining the contribution of the material to the 
Casimir pressure, the remaining contribution is the one 
associated to the EM field’s initial conditions. 

From Eg. (133)) ■ it is worth noting that unlike the other 
contributions this one does not simplify when consider¬ 
ing homogeneous and isotropic bodies. In fact, it does not 
depend directly on the material properties and bound¬ 
aries as the material’s contributions do. The dependence 
on the boundaries and the materials is encrypted by the 
EM retarded Green tensor but there is no more explicit 
dependence. 

Nevertheless, for the particular case of the Lifshitz 
problem we can perform a Fourier transform in the paral¬ 
lel coordinates through Ea. dMl) and then integrate them, 
obtaining: 


If dk 

P7c(xi,/3z.M) = -^y 


^bm _ 


k^k'^ 


w. 


k J 


X coth 




rai+ioo .a2+ioo 


27rz 


27rz 


7 (S1S2 AWk) 


X lim 

X2—>-Xi 


X J dz tz](^gj(zi,2;',k||,si) 

[ dz" , (84) 


which is a simplification that can be done when the bound¬ 
aries are parallel surfaces. 

The integrals over z' and z" are over the whole axis 
(—oo,+oo). Then, as the EM retarded Green tensor is 
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given by Eqs. (isnii-dMi) for each region, the integration sep¬ 
arates into four integrals. For example, for the first integral 
we have: 




^Ret V 
-112 




A/2 


+ 1 dz'g^^^,^^{zi,z',U,si) e 


ikz 


- 1/2 

nl/2 


^Ret,Bu \ 


+ / ^Lt,Sc(^i>2'-k||,si) e* 

J-1/2 

r+°^ / 

Jll2 


(85) 


where the first and the last are the integrations having 
the source point in each plate, while the other two are the 
integrations having the source point inside the gap. 

As all the transforms of the EM retarded Green ten¬ 
sors are given in term of exponentials, the integration is 
straightforward and can be written compactly as: 


dz' Qt,{z^,z'M,s^) 


EE 

n—1 fi 


X [(-)"-!] e 


(- 1 ) 


‘^Qz [Qz + (- 1 )” ^ikz] 


o [(-) 


n—11 


- 1 ] 




—e 




1 


.0 [( 






g[(-l)"ife^-[(-l)”-|-l]g^]| ^fj. 


H 


(n) 




D 




X H., [(-)”■'] +e„,, [(-)"] 

X (86) 


It is clear that in order to obtain the result for the 
other integral / dz" C/^((((z 2 , x", —k||, S 2 ) , we only 

have to make the replacements zi —>- 22 , si S 2 , j ^ k, 
b ^ m, k|| —k|| and —>■ —kz in Ea. (|55)) . 

Therefore, studying the analytical structure (poles and 
branch cuts) of the results of both integrals will give us 
the transient evolution and also the steady regime of the 
initial conditions’ contribution to the Casimir pressure. 
The branch cuts present in the integrands can be consid¬ 
ered by pairs (see Ref. m) and the regions in the complex 


plane where the integrands in each variable are multival- 
uated can be reduced to vertical intervals of finite length 
characterized by non-positive real parts. From Eg. (1551) we 
see that the branch cut interval with null real part cor¬ 
responds to the straight line s = i Im(s) with Im(s) G 
(—k||,k||), although it does not contribute to the steady 
state, it is important to know that it can be always con¬ 
sidered as a finite-lenght interval in the complex plane. 

Then, as for the material contributions, the key point 
to determine the long-time regime of the present contri¬ 
bution is to combine the poles in each variable in such 
a way that gives no temporal dependence. 

Therefore, we must consider complementary poles. 

At first sight, the pole at the origin in each variable 
(provided by the analytical structure of the retarded Green 
tensor) satishes the requirement to give a steady term. 

Now, to see which terms present this pole in both vari¬ 
ables, we have to consider the combinations of the terms 
of 0^^{si,S2) and (siS 2 + w^). As it happened before, 
the terms TM x TM in the integrals product in Eg. (15^ 
are the ones which give the poles at the origin in each 
Laplace variable. The different combinations finally give 
terms without poles at the origin and terms with poles of 
first and second order in both variables. 

As for the material contributions, the terms without 
poles do not contribute to the residue calculation, while 
the terms with poles of second order at the origin result 
in terms which cancel out or directly vanish. 

On the other hand, as it happened before, the poles of 
first order, in the temporal domain, correspond to deriva¬ 
tives of the Heaviside functions which gives Dirac d-functi- 
ons. This is associated to the sudden beginning of the in¬ 
teraction and then these terms give no additional physical 
information. They are only mathematical consequences of 
the description introduced for the switching-on of the in¬ 
teraction. Therefore, the poles of first order do not con¬ 
tribute either. 

All in all, the pole at the origin for both variables does 
not contribute to the pressure associated to the initial 
conditions. 

However, in this case, the poles at the origin are not 
the only ones that may contribute. The spatial integration 
over z provides additional poles which are not present in 
the material’s contribution. This poles provide the (3-1-1) 
EM generalization for the modified modes which are con¬ 
sidered as an ansatz in Refs. mu and fully demonstrated 
in Ref.[T5], for the scalar case. In the present case, from 
Eg. (1551) . it is clear that integration over z results in addi¬ 
tional denominators given by qz + {—l)'"~^ikz and -I- 
{—l)"~^ikz depending if the integration is carried out over 
the vacuum gap or over the material plates respectively. 

In fact, for the terms associated to the integrations 
over the material plates, the roots provided by the de¬ 
nominator gi"^ -I- (—l)"'“^ifcz, due to the presence of £„(«) 
inside gi”\ present negative real parts, in agreement with 
the causality properties and the fact that the field dissi¬ 
pates in such regions. For the cases where these roots are 
effectively poles of the integrand (which is not necessarily 
true as we will see below), This implies that the calcula- 
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tion of the residue will give exponential decays in time, 
resulting in terms with vanishing long-time limit. 

On the other hand, it can be immediately seen that for 
the terms associated to integrations over the vacuum gap, 
the roots provided by are s = for 

n = 1 and A :2 < 0 or for n = 2 and ^2 > 0, having null real 
parts due to the free propagation of the field inside the 
gap. This shows that s = ±fwk could be poles in general, 
because the multivalued region for the integrand is given 
by the finite-lenght interval s = i Im(s) with Im(s) G 
(—k||,k||), and these poles are always located outside it 
since we always have Wk > for every k. Moreover, as 
we anticipated before, these poles are associated to the 
mentioned modified modes for the EM field. 

However, beyond the general case, for the present one 
we only have shown that s = iioJk are roots of the denom¬ 
inator of certain terms for n = 1 and ^2 < 0 or for n = 2 
and kz > 0. Moreover, for the Lifshitz problem, the given 
form of the Laplace transform of the retarded Green tensor 
makes that s = ±fwk are also roots of the numerator for 
the same cases of n = 1 and ^2 < 0 or n = 2 and kz > 0. 
Therefore, the limits of the integrands in each variable 
when the variable tends to iiwk are indeterminate limits. 
Applying L’Hopital rule we can compute the limit giving 
a finite result. This means that, for the Lifshitz problem, 
despite s = ifcuk are roots of the denominator, they are 
not poles. Therefore, there is no residue for this case and, 
moreover, there is no steady nor transient contribution for 

S = ±JWk- 

All in all, we have proved that for the Lifshitz problem 
the initial conditions’ contribution to the Casimir pressure 
vanish at the steady state {ti —>• -l-oo): 

Pic{xi, PemJ) —>■ -Pic(/3em,0 = 0- (87) 

This result means that there are no modified modes 
for this problem that contribute to the long-time regime 
of the Lifshitz problem. Physically, the result expresses 
a relaxation process where the field dissipates during the 
transient stage until reaching a steady state in the long¬ 
time regime. Moreover, for the initial conditions’ contri¬ 
bution, the process is given by the competition between 
the dissipation of the field on the material bodies present 
and the free fluctuation in the dissipationless regions (see 
next Section for further conclusions about these features). 
For the case of the Lifshitz problem, as the material re¬ 
gions are two half-spaces while the vacuum gap has finite 
length, dissipation wins over free fluctuation and there are 
no modified modes in the long-time regime. When interac¬ 
tions begin, the transient stage take place and the initial 
free field vacuum modes start to adapt to the presence of 
the material bodies (which has also a transient dynam¬ 
ics). During this process, the free fluctuations inside the 
gap attenuate due to the dissipation exerted by the ma¬ 
terial plates and no modified modes can raise, to finally 
settle at the steady state. 

Finally, after calculating each contribution of Eqs. (EH), 
((5^ and (EH), we can write the Casimir pressure in the 
steady state {ti —>• —oo) as: 


PcU^l) PSis = PBil3B.n,l), ( 88 ) 

i.e., the bath contribution at the long-time regime is the 
one that gives the total Casimir pressure for the composite 
system. 


6 Final Remarks and Outlook 

Having solved the Lifshitz problem from a well-defined ini¬ 
tial condition problem, and achieved expressions for the 
three contributions to the total pressure for all times, we 
have shown that in the long-time limit ti —>■ —oo, the 
Casimir pressure is given only by the baths’ contribution. 
We also proved that this corresponds to the only contri¬ 
bution considered in Ref. [3], using an approach based on 
the macroscopic Maxwell equations and the source the¬ 
ory. Thus, we have presented a first principles calculation 
in the context of non-equilibrium quantum open systems 
that shows that the initial conditions are erased in the 
steady state. 

Although this looks absolutely natural at first glance, 
the results in the previous works [18111^ suggest that there 
should be important modifications for slabs of finite width. 
Indeed, for plates of finite width and for a (5-plate in 1 -|- 1 
dimensions, the contributions of the initial conditions are 
present in the steady state [13] and agree with the one 
obtained from the modified modes of Refs. mn- Looking 
back to the situations analyzed in Refs. miis], and consid¬ 
ering the results of the present paper, the whole picture 
becomes clear. It follows that a non-vanishing contribu¬ 
tion of the initial conditions at the steady state is related 
to the existence of dissipationless regions of infinite size. 
The physical interpretation is the following. When there 
are regions of infinite size where the held dissipates, it 
happens that the free huctuations in the dissipationless 
regions vanish in the long time limit because the damping 
generated by the dissipative regions overcome the free fluc¬ 
tuations. Therefore, the initial conditions’ contribution to 
the pressure also vanishes. On the contrary, for situations 
where the inhnite regions are the ones without dissipation 
(while the dissipative regions are of hnite size), the free 
held huctuations are damped only in the dissipative re¬ 
gions, resulting in modihed modes at the steady situation 
(as it happens in Ref.[T8| for the appropriate cases). 

In fact, these modes are associated to the poles ±fwk 
that result from the spatial integrations over the direction 
orthogonal to the boundaries. As we mentioned before, 
the contribution of the initial conditions can be matched 
with a contribution associated to homogeneous solutions 
in ‘steady’ quantization schemes laiiziiii], but here we are 
showing that the modihed modes of that contribution are 
the result of the dynamical adaptation of the free held vac¬ 
uum modes to the material boundaries. Therefore, on one 
hand, the modihed modes become associated to the cre¬ 
ation and annihilation operators of the initial free held. On 
the other hand, from the perspective of a ‘steady’ canoni¬ 
cal quantization scheme, the Hilbert space of the modihed 
modes in the non-equilibrium situation is spanned by the 
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same creation and annihilation operators that for a free 
field. 

All in all, we have successfully proved the foundations 
of Lifshitz theory in non-equilibrium scenarios. To reach 
this goal, we have set up a first principles quantum dy¬ 
namical problem of the EM field interacting with matter 
and subjected to uncorrelated initial conditions. We then 
solved the full dynamical problem and derived Lifshitz 
theory as the steady state by taking the long-time limit. 
The physical understanding of the contributions that enter 
in the full dynamical problem allowed us to gain power¬ 
ful insights about the occurrence of the steady state in 
different systems involving the interaction between quan¬ 
tum fields and material bodies, that commonly appear in 
Casimir physics. It would be of high interest to general¬ 
ize the results of the present work to situations of slabs 
of finite width, in order to quantify the relevance of the 
contribution of the initial conditions to the Casimir pres¬ 
sure in the long time limit. Work in this direction is in 
progress. 

As a final comment, regarding the material controversy 
around the Drude and plasma models for the conduction 
electrons, we can say that the first principle microscopic 
approach could give fruitful results to bring light to the 
discussion. As we mentioned before, here we have proved 
the non-equilibrium Lifshitz foundations for insulator ma¬ 
terial plates, with permittivity functions that result from 
the bounded electrons modeled as polarization degrees 
of freedom. Including conduction electrons at the micro¬ 
scopic model requires an action for representing their dy¬ 
namics. This may result in different analytical properties 
for the retarded Green functions, giving new contributions 
to the steady state depending on the material model con¬ 
sidered and exposing in an explicit way the limitations of 
the Lifshitz formula. This is left as pending future work. 
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